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BALLISTIC AND SUPERDIFFUSIVE SCALES IN 
MACROSCOPIC EVOLUTION OF A CHAIN OF 
OSCILLATORS 


TOMASZ KOMOROWSKI AND STEFANO OLLA 

Abstract. We consider a one dimensional infinite acoustic chain 
of harmonic oscillators whose dynamics is perturbed by a random 
exchange of velocities, such that the energy and momentum of the 
chain are conserved. Consequently, the evolution of the system has 
only three conserved quantities: volume, momentum and energy. 
We show the existence of two space-time scales on which the en¬ 
ergy of the system evolves. On the hyperbolic scale (fe _1 , xe^ 1 ) the 
limits of the conserved quantities satisfy a Euler system of equa¬ 
tions, while the thermal part of the energy macroscopic profile re¬ 
mains stationary. Thermal energy starts evolving at a longer time 
scale, corresponding to the super diffusive scaling (fe~ 3 / 2 , xe -1 ) and 
follows a fractional heat equation. We also prove the diffusive scal¬ 
ing limit of the Riemann invariants - the so called normal modes, 
corresponding to the linear hyperbolic propagation. 


1. Introduction 

Consider a chain of coupled anharmonic oscillators in one dimension. 
Denote by q x and p x the position and momentum of the particle labeled 
by x G Z. The interaction between particles x and x — 1 is described by 
the potential energy V (t x ) of an anharmonic spring, where the quantity 

U • q® Qx—i (1-1) 

is called the inter-particle distance or volume strain. The (formal) 
Hamiltonian of the chain is given by 

ft(q,p) = 5Z**(q,p), ( 1 - 2 ) 
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where the energy of the oscillator x is defined by 

e*(q,p) := y + ^(tx)- (1-3) 

The respective Hamiltonian dynamics is given by the solution of the 
equations: 

pfU 

q x (t) = ^(<hP) = ? x , (1-4) 


dH 

Px(t) = -T— (q,p), xez. 

There are three formally conserved quantities (also called balanced) for 
this dynamics 

^ ^ hx; ^ ^ Pa:; ^ ^ G (1-h) 

XXX 

that correspond to the total volume, momentum and energy of the 
chain. The corresponding equilibrium Gibbs measures v\ are parametrized 
by A = \/3,p, r], with the respective components equal to the inverse 
temperature, velocity and tension. They are product measures given 
explicitly by formulas 


dv. a = JJ exp |-/3 ^e x - pp x - rv x ^j - <7 (A) | 


dx x dp x , 


( 1 . 6 ) 


where 


and 


5( A) := 13 j + log [V2^Z(t,I3) 


Z(t,/3)-.= / exp {—/3(V(r) — rr)}dr. 

JR 


(1.7) 

( 1 . 8 ) 


The length and internal energy in equilibrium can be expressed as 
functions of A by 


r(r,P) = (t x ) Ux = j3 1 d T Q(X), p = /3 1 d p Q{\) (1.9) 


n(r, p) = (t x ) ux - V — = + (V(x x )) Ux = -dpQ{\) + rr. 

Therefore the tension r(r,u ) and inverse temperature P(r,u ) can be 
determined from the relations: 


/? = S„S(u,r), rj3 = —d r S(u,r), 


(1.10) 
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where S(u,r ) is the thermodynamic entropy defined by the Legendre 
transform 

S(u,r) = inf [—/3rr + /3u + Q(/3, 0, r)]. (1.11) 

It is expected that, after hyperbolic rescaling of space and time (fe _1 , xe~ 
the empirical distribution of the balanced quantities 

»*(*) : = [*x(t),Vx(t),e x (t)], (1.12) 

defined for a smooth function J with compact support by the formula: 

X 

converges, as e — > 0, to the solution w T (t,y ) = [r(t,y),p(t,y),e(t,y)] 
of the compressible Euler system of equations: 

d t r = d y p , 

d t p = d y r(r, u), (1.13) 

d t e = d y [pr(r, w)], 

where u(t,y ) := e(t,y) — p 2 (t,y)/2 is the local internal energy. For a 
finite macroscopic volume, this limit has been proven using the relative 
entropy method, provided the microscopic dynamics is perturbed by a 
random exchange of velocities between particles, see [9] and [3j, in the 
regime when the system of Euler equations admits a smooth solution. 
After a sufficiently long time, the solution of fjl. 13j) should converge, in 
an appropriately weak sense, to some mechanical equilibrium described 
by: 

p(x)=p 0 , r(r(x),u(x)) = r 0 , (1.14) 

where p 0 and r 0 are some constants. 

To characterize all possible stationary solutions of fll,14p can possi¬ 
bly be a daunting task. Most likely they are generically very irregular. 
But it is quite obvious that if we start with smooth initial conditions 
that satisfy (11.141) . the respective solutions of (II. 13j) remain station¬ 
ary. Also by the same relative entropy method as used in [9] and [3], it 
follows that starting with such initial profiles, the corresponding empir¬ 
ical distributions of the balanced quantities converge, in the hyperbolic 
time scale, to the same initial stationary solution, i.e. they do not 
evolve in time. On the other hand, we do know that the system will 
eventually converge to a global equilibrium, so this implies that there 
exists a longer time scale, on which these profiles (stationary at the 
hyperbolic scale) will evolve. 

There is a strong argument, stemming from both a numerical evi¬ 
dence and quite convincing heuristics, suggesting the divergence of the 
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Green-Kubo formula, defining the thermal diffusivity, for a generic un¬ 
pinned one dimensional system, see 0 0. Therefore, we expect that 
the aforementioned larger time scale (at which the evolution of these 
profiles is observed) is super diffusive. Furthermore, an argument by H. 
Spohn [TO], based on fluctuating hydrodynamics and mode coupling, 
also suggests a superdiffusive evolution of the heat mode fluctuation 
held, when the system is in equilibrium. Consequently, one can conjec¬ 
ture the following scenario: after the space-time rescaling (e“"f, e _1 x), 
the temperature T(t,y ) = /3“ 1 (t, y) evolves following some fractional 
(possibly non-linear) heat equation. The choice of the exponent a may 
depend on the particular values of the tension and of the interaction 
potential V(r). 

In the present work we prove rigorously that this picture of two 
time scale evolution holds for the harmonic chain with a momentum 
exchange noise that conserves the total momentum and energy. We 
consider the quadratic Hamiltonian: 

P) — ^ ] 2^ x ~ q y 1 a x—x'(s\x — Qx') (1-15) 

X X,x' 

where harmonic coupling coefficients a x - x >, between atoms labelled by 
x and x 1 , are assumed to decay sufficiently fast. The detailed hypothe¬ 
ses made about the potential are contained in Section l2Tl Furthermore 
we perturb the Hamiltonian dynamics, given by 01.41) . by a random ex¬ 
change of momentum between the nearest neighbor atoms, in such a 
way that the total energy and momentum of the chain are conserved, 
see formula (12.7ft below. We assume that the initial configuration is 
distributed according to a probability measure p e such that its total 
energy grows like e -1 (i.e. macroscopically, it is of order 1). Addi¬ 
tionally, we suppose that the family of measures (/4 e ) e>0 possesses a 
macroscopic profile w T (y ) = [r(y),p(y),e(y)\, with r,p,e belonging to 
C£°(R) - the space of all smooth and compactly supported functions. 
The above means that 

lh nl e '52( tVx )^ J ( ex ) = [ w T (y)J(y)dy , (1.16) 

for any test function J € G^°(r). 

We decompose the initial configuration into a phononic (low frequen¬ 
cies) and thermal (high frequences) terms, i.e. the initial configuration 
(t, p) is assumed to be of the form 

c, = c, + c and + p", x € Z, 


(1.17) 
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where the variance of (x", p") around the macroscopic profile (r(ex),p(ex)) x 
vanishes, with e —» 0+ (see (I2.28P and (12.291) below). This obviously 
implies that the macroscopic prohles corresponding to (t(,, p(.) equal to 
0. Concerning the distribution of this configuration we suppose that 
its energy spectrum, defined in (I2.23jl below, is square integrable in 
the sense of condition (12.26ft . We call the respective energy prohles 
corresponding to (t^, pj.) and (t", p") as the thermal and phononic (or 
mechanical) ones, see Definitions 12.21 and 12.31 below. In Section [9] we 
give several examples of families of initial configurations satisfying the 
above hypotheses. Among them are inhomogeneous Gibbs measures 
with the temperature, momentum and tension prohles changing on the 
macroscopic scale, see Section 19.2.31 These examples include also the 
families of measures corresponding to the local equilibrium states for 
the dynamics described by (11.21) . see Section 19.2.51 

In our hrst result, see Theorem 13.II and Corollary 13.21 below, we show 
that the empirical distribution of It ) x (t) converges, at the hyperbolic 
time-space scale, to a solution satisfy the linear version of (I1.13P with 
r(r,u ) := Tir, where the parameter ri, called the speed of sound , is 
given by formula (12.51) . 

Notice that since we are working in infinite volume, relative entropy 
methods cannot be applied to obtain such hydrodynamic limit, as done 
in [3] and [9]. Instead, we obtain the limit result, using only the tech¬ 
niques based on the L 2 bounds on the initial configurations. After a 
direct computation it becomes clear that the function 

T(y) ■■= e(t, y ) - e ph (t, y ), (1.18) 

is stationary in time. The macroscopic phononic energy is given by 

e P h(t, y) ■= P V ^ + Tir 2 (t, y ), (1.19) 

Also, the entropy in this case equals 

S(u,r) = C log(w — Tir 2 ) 

for some constant C > 0 depending on coefficients ( a x ). In fact, as it 
turns out, see Theorem 13.11 below. T(y) and e p h (t,y), given by (11.18P 
and (I1.19P above, are the respective limits of the thermal and phononic 
energy profiles. 

Concerning the behavior of the thermal energy profile at the longer, 
super diffusive time scale it has been shown in [5] that after the space- 
time rescaling (e~ 3//2 f, e~ 1 x), the thermal energy profile converges to 
the solution T(t,y) of a fractional heat equation 

d t T(t, y) = - c\A y \ 3/i T(t , y), with T(0, y) = T(y) (1.20) 
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and the coefficient c given by formula (I5.50j) below. This result is 
generalized in the present paper to configurations of the form (II. 171b 

Finally, in Theorem 13.41 we consider the empirical distributions of 
the microscopic estimators of the Riemann invariants of the linear wave 
equation system that describe the evolution of the macroscopic profiles 
of (t x (t), pz(£)) - the so called normal modes. We show that they evolve 
at the diffusive space-time scale (fe~ 2 ,r/e _1 ). A similar behavior is 
conjectured for some an-harmonic chains, e.g. those corresponding to 
the /3-FPU potential at zero tension (cf. [ID]). 

Concerning the organization of the paper in Section [2] we rigorously 
introduce the basic notions that appear throughout the article. The 
main results are formulated in Section [3j In Sections [4] - [8] we present 
their respective proofs. Finally in Section [9] we show examples of dis¬ 
tributions of the initial data that are both of phononic and thermal 
types. 


2. The stochastic dynamics 

2.1. Hamiltonian dynamics with a noise. Concerning the interac¬ 
tions appearing in (jl. 15j) we consider only the unpinned case, therefore 
we let 

Oq • ^ ^ (2.1) 

xj^O 

Define the Fourier transform of (a x ) x&z by 

0(4) :=£>* exp {— 2nikx} , k G T, (2.2) 

X 

where T is the unit torus, that is the interval [—1/2,1/2] with the 
identified endpoints. Then, from (12.ip . we have d(0) = 0. 

Furthermore, we assume that: 

al) coefficients (a x ) xez are real valued, symmetric and there exists 
a constant C > 0 such that 

\a x \ < Ce~^ c for all i6Z. 

a2) stability: a(k) > 0, k ^ 0, 
a3) the chain is acoustic : d"(0) > 0, 

The above conditions imply that a{k) = a(—k), k E T. In addition, 
a G (^“(t) and it can be written as 

a(k) = ATiS 2 (k)tp 2 (s 2 (k)). 


(2.3) 
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The function p : [0, +oo) —* (0, +oo) is smooth and </?(0) = 1. For the 
abbreviation sake, we shall use the notation 

s(k ) := sin(7r/c) and c(k) := cos(7 rk), k £ T. (2.4) 

The parameter 77 appearing in ('12.31) . called the sound speed , is defined 


by 




&"( 0 ) 

Tl 87 r 2 ‘ 

(2.5) 

Let also 




tp _ := inf p and tp + := sup <p. 

(2.6) 


Since the Hamiltonian is invariant under the translations of the po¬ 
sitions q x of the atoms, the latter are not well defined, and the configu¬ 
ration space of our system is given by ((t x ,p a; )) a ; ez , where t x should be 
thought as the distance between particle x and x — 1. The evolution is 
given by a system of stochastic differential equations 

dv x (t) = V*p x (t) dt, (2.7) 


dp x (t) = |-a* q*(i) - ^{/3 *p(t)) x }dt 

+7 1/2 {Y x+Z p x (t))dw x+Z {t), xGZ. 

z=-1,0,1 

with the parameter 7 > 0 that determines the strength of the noise in 
the system. The vector fields (Y x ) xez are given by 

Y x • (pz Px+i)^p x _i T (px+i Px— i)^pi T (px —1 PxO^pz+i - (2-8) 

Here ( w x (t ), t > 0) xgz are i.i.d. one dimensional, real valued, standard 
Brownian motions, that are non-anticipative over some filtered proba¬ 
bility space (fi, J 7 , (Jy), P). The symbol * denotes the convolution, but 
in particular, since d(0) = 0, we can make sense of (remember that q x 
are not really defined) 

a * q x := ^ ^ ol x — x i q x ' :X , 

x' 


where 


' if x' < x, 

q x ’,x := 0, if x' = x, 

k £*<*"<*'**", if T < x'. 


(2.9) 
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Furthermore, /3 X = A/3® and V*/3x°\ where 

f -4, x = 0, 

fi 0) = < -1, x = ±1, (2.10) 

0, if otherwise. 

The lattice Laplacian of a given (g x ) x&z is defined as A g x := g x+ \ + 
9x- 1 - 2^. Let also V<? x := g x+l - g x and V*g x := g x -1 - # x . A simple 
calculation shows that 

P(k) = 8 s\k) (1 + 2c \k)) = 8 s 2 (k) + 4s 2 (2fc), (2.11) 

and 

P {l \k) = 4 ie~ 2ink s(k) (1 + 2c 2 (/c)) . (2.12) 

The evolution equations are formulated on the Hilbert space £2 made 
of all real valued sequences ((t x , px)) x(S z that satisfy 

E(^ + ^) < + 00 - (2.13) 

X 

Concerning the initial data we assume that it is distributed according 
to a probability measure /i e that satisfies 

sup < +°°. (2.14) 

^ 6 ( 0 , 1 ] x 

The above means that the total macroscopic energy of the system is 
finite. Here (-) Me denotes the average with respect to /i e . 

Denote by E e the expectation with respect to the product measure 
P e = <S> P. The existence and uniqueness of a solution to (12.711 in i 2 , 

with the aforementioned initial condition can be easily concluded from 
the standard Hilbert space theory of stochastic differential equations, 
see e.g. Chapter 6 of [4], 

2.2. Energy density functional. For a configuration that satisfies 
fl2.13[) . the energy per atom functional can be defined by: 

e*(H P) := Y “ \ E ( 2 - 15 ) 

x' 

Notice that = P)- We highlight here the fact that, although 

the total energy is non-negative, in light of the assumptions al)-a3) 
made about the interaction potential, the energy per atom z x does not 
have definite sign. However we can prove, see Appendix [A] below, the 
following fact. 
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Proposition 2.1. We have 

+ - C +X( r * + P *)’ (t,p)e£ 2 , ( 2 . 16 ) 

XX X 

with C- := min{ 1/2, T\tp 2 _} and c + := max{l/2, In addition, 

c * i e ^( r ’ p)i - X e *( r ’ p ) - X mtp)i ( 2 - 17 ) 

X XX 

and c* := c_ (max {1/2, X^ >0 z 2 \a z \}) 1 . 

Thanks to (12. 14ft and (12.161) the distribution of the initial data is 
such that the macroscopic energy of the chain at time t = 0 is finite, 
i.e. that 

K 0 ■= sup < Too. (2.18) 

< 0 , 1 ] “ 

We shall also assume that there exists a macroscopic profile for the 
family of measures (/r) e>0 , i.e. functions r(y), p(y ) and e(y) belonging 
to Cg°(R) such that the respective w T (y ) = [r(y),p(y),e(y)] satisfies 

(PUSD- 


2.3. Macroscopic profiles of temperature and phononic energy. 

The macroscopic temperature profile corresponding to the family of laws 
(/i e ) eS (o.i] is defined as a function T : R —» [0, Too) such that 


f T(y)J(y)dy= hm e^J(ex) j ((^P*) 2 )^ “7 

^ x l x' 

(2.19) 

for any J G Cq°(r), where S £ p x := p^, — p(ex) and 

' £*'<*"<* for", if x' < x, 


S e q X ', x ■= 0, if x' = x, 

'E x< x"<x>fex», ifx<x', 


( 2 . 20 ) 


with S € t x := x x —r(ex). After a simple calculation one gets the following 
identity 

T(y) = e(y) - e ph (y), (2.21) 

with 

e P h (y) ■= - ( p 2 (y ) T Tir 2 (r/)) , ye R, 
the phononic (or mechanical ) macroscopic energy profile. 


( 2 . 22 ) 
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2.4. Energy spectrum. The energy spectrum of the configuration 
distributed according to /z £ is defined as 

m«W := (IpWI 2 ^ + WL{ m \\, (2.23) 

where p(k) and r (k) are the Fourier transforms of (p x .) and (t x ), respec¬ 
tively. Using (12.3p we get 

tT) e (fc) = (|p(fc)| 2 >^+Ti^ 2 (s 2 (fc))(|r(A:)| 2 ) /iE , ke T. (2.24) 
Assumption (j2.18j) is equivalent with 

K 0 = sup e / t v t {k)dk < +oo. (2.25) 

£6(0,1] J T 

Definition 2.2. The family of distributions (/a e ) e>0 is said to be of a 
thermal type if its energy spectrum ttp(fc) satisfies 

K\ := sup e 2 / rv 2 (k)dk < +oo. (2.26) 

£6(0,1] J T 

Remark. In Section [5. 1 I we show that (12.261) implies that the respec¬ 
tive macroscopic prohles {r(y),p(y)) vanish. Therefore, we conclude 
that then the macroscopic phononic energy e p h (y) = 0 (see (I2.22p l 
and, as a result, 

limeJ2 J (^)M„=f J (y) T (v) d V, JeQ). (2.27) 
e_>,u+ x Jr 

We stress that although r(y) = 0 and p(y ) = 0, in case condition 
(I2.26P holds, conhguration (r, p) need not necessarily be centered (in 
the measure /i e ). Condition (12.261) is related to the issue of variability 
of the initial data on the microscopic scale, as can be seen in examples 
presented in Section [91 

Definition 2.3. The family (/i e ) is called to be of a phononic (or 
mechanical) type if 

(t r - - r ( ex ^ 2 )^ = 0 ( 2 - 28 ) 

X 

and 

= °> ( 2 - 29 ) 

X 

where r(x) and p(x) are the macroscopic prohles of the strain and 
momentum, see (11.161) . 
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Remark. Using (12. 19ft and (12.28ft together with (12.29ft it is straight¬ 
forward to see that the temperature profile corresponding to a phononic 
type initial data vanishes, i.e. T(y) =0. 


3. Statement of the main results 

3.1. Hyperbolic scaling. Assume that the configuration (r, p) can be 
decomposed into two parts whose laws are respectively of the thermal 
and phononic types. More precisely, we assume that 

r 3! — t x + t ac, Px — p'x + Pi) xez, (3.1) 

where 

K[ := sup e 2 j(tv' e ) 2 (k)dk < +oo, (3.2) 

eG(0,l] J T 

and ro' (/c) is the energy spectrum of (t 7 ,p 7 ), i.e. 

K(k) = (Ip'WI 2 )^ + n <p 2 (s 2 (/c)) , k e T. (3.3) 

with (p'(fc), r'(/c)) the respective Fourier transforms of (p 7 ,t 7 ). In addi¬ 
tion, we assume that the configuration (t 77 , p 77 ) is of the phononic type 
in the sense of Definition 12.31 

Suppose that (t 7 (f),p 7 (f)) and (t"(f), p"(£)) describe the evolution of 
the respective initial configurations (t(,,p(.) and (t",p") under the dy¬ 
namics (12.71b Define the microscopic temperature and phononic energy 
density profiles by 

eth.zfa) := ^(t 7 (t),p'(t)) and e phiX (t) := e x (t"(f),p"(t)) (3.4) 

Let 

v T (t,y) := (r(t,y),p{t,y)) (3.5) 

be the solution of the linear wave equation 
{ d t r(t,y) = d y p(t,y), 

l dtP(t,y) = r 1 d y r{t,y), (3.6) 

[ r(0, y) = r(y), p(0,y)=p(y). 

The macroscopic phononic energy at time t, see (|2.22lh is given by 

e P h (t, y) = \ (p 2 (^ y) + nr 2 {t, y)) ■ (3.7) 

The components of 


w T {t,y) := [r(t,y),p(t,y),e ph (t,y)] 
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evolve according to the system of linear Enlcr equations: 

( d t r(t, y) = d y p(t, y), 

I d t p(t,y) = T 1 d y r{t,y), 

d t e ph (t, y) = T X d y (r(t, y)p(t, y )), 

[ w(Q,y) = w(y). 

The following result is proven in Section [6j 


Theorem 3.1. Suppose that the initial configuration (r, p) satisfies the 
assumptions formulated in the foregoing. Then, (r'(te _1 ), p'(te -1 )) is of 
thermal type for all t > 0 and 

/ +°o f t\ r+oo r 

J(t,ex)E e t t h, x ( - j dt = / dt T(y)J(t,y)dy 

(3.9) 

for any J G Co°([0, +oo) x R). Configurations (r "(te 1 ),p // (te x )) are 
of phononic type for all t > 0 and 


lim e > 

e—>0+ ^ 


A p h. X ( 


®ph ifi 


ex 


= 0. 


(3.10) 


In addition, 



J{t,ex) E e c x 



for any J G Qf’QO, +oo) x R). 


[e P h(t, y) + T(y)\ J ( t , y)dtdy, 
(3.11) 


Concerning the macroscopic evolution of the vector tv x (fi), see fll,12p . 
the above result implies the following. 

Corollary 3.2. Assume that 

^ e\] J (ex)(w x (0)) fle = I J(y)w(y) dy, J G C“(R). (3.12) 

+ x 

Then, for any J G C^°([0, +oo) x R) we have 

r+oo f t\ r+°o r 

lim e y / J(t,ex) E e tu x ( - ) dt = / J(t,y)w(t,y)dtdy , 

e _>°+ z^j Q W j Q j r 

(3.13) 

with the components of w(t,y) being the solutions of the linear Euler 
system (13.81) . 
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3.2. Super-diffusive scaling. It follows from Theorem 13.11 that the 
phononic energy evolution takes place on the hyperbolic scale that is de¬ 
scribed by the linear Euler equation. In consequence, it gets dispersed 
to infinity on time scales longer than £e _1 . On the other hand, it has 
been shown, see Theorem 3.1 of [5j, that the respective temperature 
profile T(y), see (I2.19|) . evolves on a superdiffusive scale (te~ 3 / 2 , xe _1 ). 
Therefore, we have the following result. 

Theorem 3.3. Suppose that the distribution of the initial configura¬ 
tion (t, p) satisfies condition (13.11) . Then, for any test function J e 
C“([0, +oo) x R) we have: 

/ +°c / t \ r+oo r 

J(t, ex)E e c x ( ) dt = / / T(t,y)J(t,y)dtdy. 

(3.14) 

Here T(t,y ) satisfies the fractional heat equation (15.49j) with the initial 
condition T(0,y) = T(y), given by (12.27j) . 

The proof of the above result is presented in Section d 


3.3. Evolution of the normal modes at the diffusive time scale. 

Finally we consider the Riemann invariants of the linear wave equation 
system (13.61) . They are given by 

■= P(t,y) ± y/nf{t,y). (3.15) 


The quantities defined in (j3.15j) are constant along the characteristics of 
(13. 6 p . which are given by the straight lines yiyTR = const. Therefore 

=/<*>(!, ±VF«), (3.16) 

with / - qy) determined from the initial datum. This motivates the 
introduction of the microscopic normal modes given by 


fi +) W :=P y(t) + y/n 





Xy(t), 


f y \t) '■= Py(f) - Vn 


1 


1 

2 






(3.17) 


that are the second order approximations (up to a diffusive scale) of 
f^(t, y). The particular form of f^(t) is determined by the fact that 
we are looking for quantities of the form p(t)±^/ri (1 ± c±V*) x(t) that, 
at the hyperbolic scale, approximate the Riemann invariants (j3. 15|) 
and possess limits at the diffusive time scale. The latter requirement 
determines that c± = (1/2) ( 37 /yTi ± l). 
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Remark. The normal modes fW(t) capture the fluctuations around 
the Riemann invariants of the linear hyperbolic system (1 1.13 [1 . Their 
analogues, called normal modes (j>±(t), are considered in Appendix 1 
b) of m, i n the context of anharmonic chains with FPU-potential, see 
(8.18) of ibid. 


To state our result rigorously, assume that (p e ) e - the initial laws of 
configurations (t, p) satisfy (jl . 16j) . Denote also the heat kernel 


Pt(v) ~ 


: exp 


y 2 ) 

' ADt J ’ 


(3.18) 


\J AixDt 

where D := 3y. Let fi d \t,y) := P t * fd\y), i = ± and 

f {± \y) := p{y) ± Vn r (y)- (3.19) 

With the above notation we can formulate the following result. 

Theorem 3.4. Suppose that condition (11.161) is in force. Then, the 
phonon modes fy^it) satisfy 


f 1 i 1 ? + e E J ( ex ) E ^ ±) (“) = 

for any J G Cq°(r). In addition, for any t G {—, +} 


y)dy. 


(3.20) 


e 1 T + e E J ( ex_t v / u;) E e f; 


p(0 

fx 


= / f!: d \t,y)J(y)dy, (3.21) 


The proof of this theorem is contained in Section [8j 


Remark. The results of the present paper are valid also for other 
stochastic dynamics obtained by a stochastic perturbation of the har¬ 
monic chain that conserves the volume, energy and momentum. For 
example we can take a harmonic chain perturbed by a noise of the 
’’jump” type. More precisely, let (N XtX+ i(t)) xez be i.i.d. Poisson pro¬ 
cesses with intensity 3y/2. The dynamics of the strain component 
(br(t)) xez is the same as in (12.7jl . while the momentum (p x (t)) xgz is a 
cadlag process given by 

dp x (t) = —(a * q(t)) x dt 

(3.22) 

+ [Vpo;(t— )dN XjX+ i(t) + V*p x (t— )dN x _i^ x (t )], x G Z 

where (a * q) x (t) is defined as in (12.91) . One can show, with essen¬ 
tially the same arguments as the ones used in the present paper, that 
Theorems 13.11 - l3~4l hold also for this dynamics. 
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4. Asymptotics of the phononic ensemble 


In this section we assume that the laws (/x e ) e >o of the initial con¬ 
figuration (t, p) is of the phononic type, i.e. they satisfy (12.28P and 

tra. 


4.1. Approximation of the macroscopic phononic energy. The 

solution of (13 .6j) is given by 


r(t,y ) 


= r(y + Vnt) + r (y — 0T t) + 


Ip (y + Vnt) - p(y - Vnt)\ 

\/a 


(4.1) 

p(t, y)=p(y + Vrt) +p(y- y/nt) + Vn[r{y + y/nt) -r(y- i/rlt )], 


where r(y),p(y) are the initial data in (1.3.61) . Using the above formulas 
and (13.7p we can write e p h (t, y) in terms of the profiles r(y) and p(y). 
We can easily infer the following. 


Lemma 4.1. Suppose that r,p G Co°(r). Then, 


lim sup 

e -> 0 + t> o 


e£e ph (t,ex) - / e ph (t, 

Jr 


y)dy 


= 0. 


(4.2) 


4.2. Evolution of the mean ensemble. Define the mean configura¬ 
tion 


where 


4 e) 0) : = e*(t £ (t),p e (t)) = 


put) 


(f))iez, 


= Eepz(t)- 

(4.3) 

n as 


‘_j a x—x'Ae,x,x'(t)i 

(4.4) 


x' 

where q e ,x,x'(t) is the respective mean of q x y(t) defined by (12.9p . 

Let 6 1 (t,k) := [i(t, k), p(i, k)], k G T, where i(t,k),p(t,k) are the 
Fourier transforms of the respective components of the configuration 
(r x (t), p x (£))). Suppose that 5 > 1. Define 

x e (t, k) := eE e r (Jj, elJj and jj e (t, k) := eE e p ek^j . (4.5) 

Conservation of energy and condition (12.18P imply that 

E* := sup f (\t e (t,k)\ 2 + \$ e (t,k)\ 2 )dk <+oo. (4.6) 

t>0,ee(0,l] J e _1 T 
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From (12.7ft we obtain 


d D £ (t, k) = k), 


dt 


where 

and 


= [i ,(t,k),p,{t,k)} 


A,{k) := 


0 


(4.7) 

(4.8) 

(4.9) 


—r(efc)a* —b 

with a* - the complex conjugate of a and 

f(k) := Ti<^ 2 (sin 2 (7r/c)), a := -(1 — exp {— 2mek}), 

b := — 11 — exp {—2niek} [ 2 [2 + cos(27t ek)\. 

Observe that for a given M > 0 we have 

lim sup \A e (k) — A 0 (k)\ = 0 

e ” >0 + \k\<M 


(4.10) 

(4.11) 


and 


Mk) ■= 


0 


2irik 


2niTik 0 

Since both eigenvalues of Ao(k) are imaginary we have 

D * := sup || exp {A 0 (k)t} || < +oo. 

(t,fc)e r 2 


(4.12) 


Here || ■ || is the matrix norm defined as the norm of the corresponding 
linear operator on a euclidean space. In fact, an analogue of (14.121) 
holds in the case of the linear dynamics governed by (14. 7 j) . 


Lemma 4.2. We have 


C * := sup || exp {A e (k)t} || < +oo. 

(t,fc)GR 2 ,ee(0,l] 


(4.13) 


Proof. The eigenvalues of A e (k) equal 

A± = 16 ± \/b 2 - 4f(e/c)|a| 2 | . 

Since both b and f (ek) (see (j4.10|) ) are real and non-negative it is clear 
that ReA-t < 0 and the conclusion of the lemma follows. □ 
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4.3. Asymptotics of the mean ensemble. Let x>^\t,k) be the so¬ 
lution of the following equation 


(J, k) = Aj A 0 (kpV(t, k), 6<°>(o, k ) = s.(o, k). (4.14) 

Let also v(t,£) be the Fourier transform of v(t,y ) the solution of the 
linear wave equation system given in (13.51) . It satisfies 

v(t, £) = exp {A 0 (£)t} 7(0, £). (4-15) 

Conditions (I2.28P and (I2.29|) imply 


lim / 

e-i> 0 + J e -lj 


o e (0,£) - v{0,£) 


d£ = 0. 


Lemma 4.3. Suppose that S 6 [1,2) and M > 0. Then, 


'>M 


lim / |o e (f, k) — V^ 0) (t, k)\ 2 dk = 0, t > 0. 
^o+ J__ M 

Proof. Note that A e (k) = A 0 (k) + eB e (k ) (see (14.91) ). where 
K(M) := sup sup |||| < +oo. 

£6(0,1] \k\<M 


(4.16) 


(4.17) 


(4.18) 


From (14.71) we can write 

^B e (t, k) = -J^A 0 (k)v e (t, k) + e 2 ~ s B e (k)£i e (t, k ), (4.19) 

Therefore, by Duhamel’s formula, we conclude 

b e (t, k ) — 6®(t, k) = e 2 ~ 5 J exp |^—pA 0 (fc)| B e (k)t) e (s, k)dk. 

(4.20) 

The lemma then follows from (14.61) . (j4.18j) and Lemma [4.21 □ 


Lemma 4.4. For any t > 0 and 6 G [1, 2) we have 

2 


lim / 

e-A)+ J e -i T 


* e (t,£)-v( — v £ 


t 


and, in consequence, 

lim e 

6—^0+ 


M ft) 


d£ = 0 


= 0. 


(4.21) 


(4.22) 
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Proof. Formula (14.22ft follows from (14.21ft so we only focus on the proof 
of the latter. We have 


£) = exp 



Oe(0,^). 


Therefore,from the above and (14.151) we can write 


limsup / 

£—>■0+ J e~ 1 T 




d£ 


< 2 lim sup / 

e->0+ Je~ 1 T 

+2 lim sup / 
£->■ 0 + J e -lT 


expMeW-jzr ? O e (0,£)-f>(0, 


d£ 


(4.23) 


exp A e {£)-£-£ V v(0,£) - exp A 0 (£)- s — [• v(0,£) 


d£. 


Using Lemma EOl and (j4.16j) we conclude that the first term on the right 
hand side vanishes. To estimate the second term divide the domain of 
integration into the regions \£\ < M and its complement. On the first 
region we use Lemma 14.31 to conclude that the respective limit vanishes. 
Therefore, we can write that, modulo multiplication by factor 2, the 
second term on the right hand side of (14.23ft is equal to 


limsup / 

£->0+ J M<\t\<l/(2e) 


exp A e (£)-j—^ ^ u(0, £) - exp > v(0, £) 


d£ 


< 2 lim sup / 

£->0+ J M< |f|<l/(2£) 


exp A e (£)-j—^ U(0,£) 


+ 


exp Aq(£)~— U(0,£) 


< 4 (C 1 * + DA) limsup / |f)(0, f?)| 2 d£, 

£ ^ 0 + J M<\£\<l/(2e) 


the last estimate following from Lemma [4.21 and estimate (14.12[) . We 
can adjust M > 0 to become sufficiently large so that the utmost right 
hand side of (14.24)1 can be arbitrarily small. The conclusion of the 
lemma therefore follows. □ 


Suppose also that (■ r(t,y),p(t,y )) is the solution of the linear wave 
equation (13.bp and e p h (t,y) is the corresponding macroscopic phononic 
energy density profile, defined by (13.71) . As a direct corollary from 
Lemma 14.41 we conclude the following. 



(4.24) 
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Corollary 4.5. For any t > 0 we have 

' > te,x I 1 
6^0+ 7 ^ 


r(t, ex) 


= 0, 


(4.25) 


lim e > 

£->0+ 


1 2 


- - p(t, ex) 


= 0. 


Concerning the behavior of energy functional tx ( t) corresponding 
to the average phononic ensemble (see (14.4p ) we get the following. 


Corollary 4.6. For any t > 0 and 6 e [1, 2) we have 

t 


lim e > 

e->0+ 


>(d I L 


6ph 


-5-1 


, ex 


= 0, 


(4.26) 


where e p h (t, y) is given by (13.71) . In addition, for any 6 G (1, 2) we have 
Um =0, J e Cg°(R). (4.27) 


Proof. In light of Lemma 14.41 only equality (14.27ft requires a proof. As¬ 
sume that 5 G (1,2). Note that the left hand side of (14.27ft equals 
lim e ^ 0+ [j} £ \t) + ji c \t)}, where 


Ji%) :=iEP 

1 /-V(2e) 


ox \ e s 


J(ex) 


(4.28) 


— 1/(2e) 


J e (£) 


»l/(2e) 


—l/(2e) 


—fc p £ -^,/c dfc 


and 


^2 C, W : = y 


x,x' 


-£ 


T 


•l/(2e) 


' —l/(2e) 


J e (£) 
(4.29) 


" 1 /( 26 ) 


x 

where 

and 


- 1 /( 26 ) 


h z (-e(k + I))h z (ek)x e 


-5-1 : 


k ) t e I dk i 


, ... exp{-27r//c2;} - 1 , 

h z (k) ■= _ r , z6Z,l:6T 


exp {— 2irik} — 1 
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Recall that D € (t, k ) and vf\t, k ) are given by (14.71) and (14.141) respec¬ 
tively. Using Lemma 14.41 we conclude that 


}^- J ^W _ = 0 and ^ ~ ^ = °> 


rW 


r{e). 


rW/ 


where 


J, w (i) := / I p(f r; ,-f-k)p(f rT .k)dk (4.30) 


and 




r<5-l ' 


r | —^-r, — £ — k] r (k ) d£dk. 


(4.31) 


Recall that 


= ^exp v(fc),e 2 ^ ^ 

where e^ = [0,1] and (•, -) c d is the scalar product on C d . Therefore 


j} e \t) = / J(£)d£ I (f)(-fc-0,exp<j A£(-k 


c& — 1 


e 2 


h(fc),exp |t!q e 2 ^ dk. 


By an elementary application of the Riemann-Lebesgue theorem we get 


£-> 0 + 

This ends the proof of (14.27p 


lim Jj £ \t) = °, j = 1,2. 


□ 


4.4. Evolution of the conserved quantities at the hyperbolic 
scale. Our goal in this section is to prove the following result. 


Theorem 4.7. Suppose that (12.281) and (12.291) are in force. Then, 

1 2 


lim e > 

£—>0+ 


lim e > 

£-5-0+ ^ 


r(t , ex) 


pit, ex) 


= 0, 


= 0 


(4.32) 


and 


lim e > 

£^0+ ^ 


E f C-7 


Cph {t, 


ex 


= 0, t > 0. (4.33) 
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Proof. Equality (I4.33H is a consequence of (14.32ft . We only need to sub¬ 
stantiate the latter. Let (5x £ (t), 5p e (i)) := (Sx £jX (t), Sp £tX (t)) xGZ , where 

5x e>x (t) := r £jX (t) - v e , x (t), Sp £ , x (t) := p £ , x (t) ~ Pe,x(t), (4.34) 

with t e , x (t) and p £>:r (f) given by (I4.3j) . We have (see (14.4h ) 

E e e x (t e (t), p e (f)) = E £ t x (5x £ (t), Sp £ (t)) + 4 e) (t) ■ (4.35) 

By conservation of the total energy we obtain 


= e 


^2e £ x x ( 5x £ ( f -) ,5p e ( f - 




:(«) 


> e ^ E e e x (<St e , 5p e 


e^2e ph (t,ex) - 


,(d 


(4.36) 
e P h (t,ex) 


Letting e —> 0+ and using assumptions (12.28(1 and (12.29(1 we conclude 
that the limit of the energy functional appearing on the the utmost left 
hand side equals 

/ e ph (0 ,y)dy = / e ph (t,y)dy, t e R, (4.37) 

Jr Jr 

with e p h (t,y) given by (13.71) . Using the above together with (j4.26jl we 
conclude that 

e i - MM - ) ) + / e P h (t,y)dy. 


/ e ph (0, y)dy > lim sup e V' E e e x (<5r, 

Jr €-*>0+ ^ \ 

From (I4.37|) we infer therefore that 

lim sup e ^2e £ c x (<5t e ( - MP t 


I R 


£ -S>0+ 


= 0. 


(4.38) 


From (12.161) applied to the configuration (§x £ (te 1 ),<5p e (te x )) we 
conclude that there exists c* > 0 such that 


c_e 


£ : 




< eJ^E e e x f Sx £ ( - J , tip 


+ [Sp 

t 


(4.39) 


Equalities (14.32)1 then follow straightforwardly from (14.39)1 and Corol¬ 
lary 14.51 □ 
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4.5. Evolution of energy density at the superdiffusive scale. In 

this section we show that the energy density disperses to infinity at 
time scale t/e 5 , where 5 G (1,2). 

Theorem 4.8. Suppose that (12.281) and (j2.29f) are in force, and 5 G 
(1,2). Then, 

e J ( ex ) E e*x =0, (4.40) 

X ' ' 

for any J G Co°(r), t > 0. 

Proof. In analogy to (I4.36P we write 

e S = e £ ( r (75 ) > P f 4 


> e^E e e x ( 5r, 


+e e P h 


,^P £ 


(4.41) 


,< 5-1 


, ex 




>(d 


Cph 


,< 5-1 


, ex 


Therefore taking the limit, as e —> 0+, in (14.41)1 we obtain (using 
and (14.26)1 ) 

J e ph (0 ,y)dy > lim sup e ^ E e e x (^j ,Sp e (J^ 

+ limsup J e p h dy. 

By (14.37P we conclude 


(4.42) 


lim sup e E e e x I 5t, 


e->0+ 




= 0. 


(4.43) 


Since, for any J G C^°(r) we have 

e Z J ( a )E* f r 


(4.44) 


e^2j{ex)E e t x [St J-j] (-j) j +tJ2 J ( 


[ex)^ | — 


The conclusion of the theorem follows directly from (I4.43P and (14.27P . 

□ 
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5. Evolution of a thermal ensemble 

In this section we investigate the limit of dynamics of ensemble fami¬ 
lies whose initial laws (y e )e>o satisfy condition (I2.26p . Then, as it turns 
out, the respective macroscopic profile of (t, p) is trivial, i.e. r(y ) = 0 
and p(y) = 0. In fact, we show that condition (12.26[i persists in time 
and the temperature profile remains stationary at the hyperbolic time 
scale. 


5.1. Wave function. Define the wave function of a given configura¬ 
tion as 

Vtr • (*V * * 1 ) 3 ; T ipx ^ ^ ^x—x'Po,x' + ^Pxi XEZ, (5.1) 

x' 

where q 0 y is given by (12.9ft and ( Cj x ) are the Fourier coefficients of the 
dispersion relation defined as, see (12.2|h 

Lu(k) := d 1//2 (fc) = 2 ^/ti \s(k)\<p(s 2 (k)). (5.2) 

Using (15.2p we can rewrite (15.ip in the form 

Vh = (^ (1) * + ipx-, x G Z, 

with (d)^) are the Fourier coefficients of 

uj^\k) = —ie mk sign(k)y/r[(p(s 2 (k)), k G T. (5.3) 

The Fourier transform of the wave function ( ijj x ) can be written 
as 

, 4>(k) = u^\k)i(k) + ip(k). (5.4) 

Therefore, 

(llWI 2 )n, = {\u m (k)i(k)\ 2 + |p(*)|')„. + 2(Im (i(k)u m (k)p’(k))) ll .. 

(5.5) 

On the other hand, since i(—k) = i*(k) and p(— k) = p *{k) we obtain 

(l^(-m 2 )c. = (|w (1) Wt(*:)| 2 + |p(fc)| 2 )p, - 2(Im (i(k)J 1) (k)p , (k))) lu . 

(5.6) 

From (15.51) and (15.6p we obtain 

[rv 2 e (k)dk= [(\m\%dk. (5.7) 

J T J T 

Condition (12.261) therefore is equivalent with 

sup e 2 / (|^(/c)| 2 ) 2 e dfc < + 00 . 
e£(0,l] Jt 


(5.8) 
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Remark. Note that the macroscopic profile (r(y),p(y)) correspond¬ 
ing to a given configuration (r, p) satisfying (12.26j) necessarily vanishes, 
i.e. 

lim c y^(t x ) Me J(ex) =0 and lim e^(p x )^J{ex) = 0 (5.9) 


e-S>0+ 


for any J G Cq^R). To show (15.9j) it suffices only to prove that 

= °- 

e—^0+ z —' 
x 

Indeed, from conditions (12.281) and (I2.29p . we obtain 

ym = / i li(y)J(y)dy, 


(5.10) 


(5.11) 


where 

and 


^(y) = y/nH{r){y) + ip(y) 


n{f){y) := —p.v. 


7T 


[MW 

I R y-y’ 


is the Hilbert transform of a given function / G C'^°(r). The integral 
on the right hand side is understood in the sense of the principal value. 
In particular, = 0 implies that both p(y) = 0 and r(y) = 0. 

To show (15.101) observe that by the Plancherel identity we can write 
that the absolute value of the expression under the limit in (15.10[) equals 


where 


${k))i^J t {k)dk 


_^ 1 ^ / A; 

J e (k) := J{ ex ) ex P {—27 t ikx} ~ - J ( — 


(5.12) 


and J(k) is the Fourier transform of J(x). Expression in (j5. 12j) is 
therefore estimated as follows 

3/4 


{m)J 


dk 


< 


(\m\Xdk 


\_J T 


1/4 


J 


4/3 


dk 


where the estimate follows by Holder inequality. Using the change of 
variables k' := k/e in the second integral on the right hand side we 
conclude that it is bounded by eiF^ 4 || J|| l 4 /3( R ) for e G (0,1], which 
proves (15.101) . 
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5.2. Evolution of the energy functional at the hyperbolic time 
scale. Define by tv e (t, k) the energy spectrum (12.231) corresponding to 
the configuration (r(t),p(t)). The following result asserts that there is 
no evolution of the macroscopic temperature profile at the hyperbolic 
scale. Moreover, if the initial energy distribution is of the thermal type 
it remains so at this time scale. 


Theorem 5.1. For any J G C^°([0,+oo) x R) we have 


r+° o /f\ r+oo /• 


J(t,y)T(y)dtdy. 

(5.13) 


Moreover, at any time t > 0 the energy spectrum satisfies 

t 


K x (t) := sup e 2 / to 2 [ -,k ) dk < +oo 
£6(0,1] J T 


(5.14) 


and 


lim e ^ E e t ^ - j J(ex) = 0 and lim e ^ E e p ^ - j J{ex) = 0 

(5.15) 


for any J G Cq°(r). 

The proof of this result is given in Section 15.51 


5.3. Wigner transform. Define ( t ) :—if x (te x ) and k ) its 

Fourier transform (belonging to L 2 (t)). By W e (t) we denote the (av¬ 
eraged) Wigner transform of (t ), see (5.9) of [5j, given by 


(W e (t), J) [ W € (t,p, k)J*(p, k)dpdk, 

J RXT 


(5.16) 


where Wjt. p, k) the Fourier - Wigner transform of fi>^ e \t) given by 


w,(t,p,k) := |e, [(■/><•>)* (t,k - Jf) #'> (t,k + |)] (5.17) 


and J belongs to S- the set of functions on R x T that are of C°° class 
and such that for any integers /, m, n we have 


sup (l + y 2 )"|c^rA!/A)l<+<x>- 

ye R, ke T 


Let A be the completion of S under the norm 

\\A\a-= / sup \J(p, k)\dp. 

Jr k 


(5.18) 
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In what follows we shall also consider the Fourier - Wigner anti¬ 
transform of ^ e \t ) given by 

e) (t, -k + 

From the Cauchy-Schwartz inequality we get 

Thanks to energy conservation property of the dynamics and Cauchy- 
Schwartz inequality we get 

sup sup(||F e (f)|U + ||W e (t)|U') < 2/1 o, (5.20) 

eG(0,l] t> 0 

where K 0 is the constant appearing in condition (I2.25jl . As a direct 
consequence of the above estimate we infer that the family (W / e (-)) eg ^ 0 ,j 
is *—weakly sequentially compact as e —> 0+ in any L°°([0, T]; *4'), 
where T > 0, i.e. for any e n —> 0+ one can choose a subsequence 
e n / such that (W e , (•)) , that is *-weakly converging. In fact, using 
hypothesis (12.261) one can prove that, see Proposition 9.1 of [5], the 
following estimate holds. 

Proposition 5.2. For any M > 0 there exists C\ > 0 such that 


IK 1 


t,k + 


ep 


(5.19) 


Y,(t,p,k) := -E, 


•t/il 


l|lk(*.P. Ollwcr)+IIU(t,P,Olll»(T) (5-21) 

< (||P(0,P, •)&<•,) + llk(0,p,-)lli=(„) e c '“, v £ e (0,1], |p| < M, t > 0. 

The right hand side of (I5.21[) remains bounded for e G (0,1] thanks 

to (J52D- 

In fact asymptotically, as e G 0+, the function | (t) | 2 coincides 
with the energy density t x (te -1 ), which can be concluded from the 
following result, see Proposition 5.3 of [5]. 

Proposition 5.3. Suppose that condition (15.8j) holds. Then, 


lim e > J(ex) E, 

e—>0+ ^ 




o, t > o, J e C 0 °°(R). 

(5.22) 
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From the above result we conclude that for a real valued function 
J(y, k ) = J(y) we can write 



X 


= lim e > J(ex)E e t x 

f->0+ ^ 

X 



(5.23) 


5.4. Homogenization of the Wigner transform in the mode 
frequency domain. It turns out that in the limit, as e A 0+, the 
Wigner transform W £ (t,x,k) becomes independent of the k variable 
for any t > 0. To avoid boundary layer considerations at the initial 
time t = 0 we formulate this property for the Laplace-Fourier transform 
of the respective Wigner function. More precisely, let 

r+oo _ 

w Ct ±(X,p,k) / e~ xt W e {t,p,±k)dt 1 (5-24) 

J o 


W,±(A,p, k) 



e xt U e ^±(t,p,k)dt, 


where 


Ue,+ (t,P, k ) 


Ue,-{t,P, k ) 


Y e (t,p, k) +Y*(t, —p, k) 
Y e (t,p, k) — Y*(t,—p, k) 


Thanks to Proposition 15.21 the Laplace transform is defined for any 
A > 0 and for any M > 0, compact interval I C (0, +oo). In addition, 
we have 


Ci := sup sup 

e€(0,l] \ei,\p\< M 


^||w e , t (A,p)|| L 2( T ) 


Let 


IKi(A,p)IUa (T) 

tG{ —,+} 


< Too. 
(5.25) 


U 'e ±) (A,P) ;= (We,+ (. A),C±) L 2 (t) = (w e -(\), C±)l 2 (t) , 

where 

c + (k) := ^s 4 (fc), C-(k) := 2s 2 (2 k). 

o 

The following result holds, see Theorem 10.2 of (5j. 


(5.26) 
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Theorem 5.4. Suppose that the initial laws satisfy (12.26ft . Then, for 
any M > 0 and a compact interval I C (A 0 , +oo) we have 

lim sup / \w e:+ (X,p,k) — w < f k \X,p)\dk = 0, (5.27) 

e_>0+ \ei,\p\< M Jt 

and 

lim sup / |w e ,±(A,p, k)\dk — 0. (5.28) 

e_s>0 + \ei,\p\<M J t 

5.5. Proof of Theorem 15. 11 Applying (15.21ft at p = 0 we conclude 
that 

sup e 2 f (e e \'ip ( ' e \t,k)\ 2 ') dk < +oo. (5.29) 

Condition (15.14ft then follows directly from (|5.7[i . In addition equalities 
(15.15[) can be inferred from (I5.29[) and (15.9ft . 

Concerning the proof of the convergence of the energy functional in 
(15.13[) . recall that we already know that (W e (-)) is *—weakly sequen¬ 
tially compact in (L 1 ([0, +oo), .4))*. Therefore for any e n —> 0, as 
n —> Too, we can choose a subsequence, denoted in the same way, such 
that *—weakly converging to some W 6 (L 1 ([0, +oo),^4))*. 

In light of (I5.20p we have 

supHWe^lU < K 0 , (5.30) 

t> o 


with Kq the same as in (12.25ft . Therefore, the respective Laplace- 
Fourier transform w e (X,p, k) can be defined for any A > 0. To identify 
W it suffices to identify the *-weak limit in A! of the Laplace transforms 
■uy n (A), as n —* +oo. Thanks to Theorem 15.41 any limit w(\,p,k) = 
w(X,p) obtained in this way will be constant in k. In fact 


w(X,p) 


lim 

n—>-+oo 


^6?( A >P)- 


We claim that in fact for any J G Cq(R) 


(w(X),J)= [ J(y)T(y)dy , (5.31) 

Jr 

which implies that in fact W G L°°([0, +oo), A!) and 

W(t,y, k) = T(y), t> 0. 


In light of Proposition 15.31 this allows us to claim (15.13ft . The only thing 
yet to be shown is therefore (15.31ft . 
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Proof of (15.31 j) . Letting = ^±(0) and uf± = U e} ±(0) we obtain 
that for any A > 0 (see (10.5) of [5]j 


Xw e>+ — = —i5 e ojw e:+ — i^R'pu e _ + — C (mi £ i + — u t)+ ) + erf\ 

(5.32) 


p r (0) 7 

A u e + — U e 1 — — u e _ H —L 
e ’ e 


U e,+ ~ ^ K,+ + W e,-) 


+ erf\ 



2u> _ 

' U e,+ 




w e>+ ) + erf ] 


\w e - — wfl = iS e ujw e _ + i~/R'pu e ^ + —C (w e ,_ — m £ ,+) + erf\ 


7 


where 


S e u(k,p) ■= - 
e 


ep 


uj [k + — ) — u [k — — 


ep 


u{k,p) \= 


2 L 


P 


u [k H— I + c o [k — 


p 


(5.33) 


Cw{k) := 2 / R(k, k')w(k')dk' — 2R(k)w(k), w G L 1 (t), 


and 


R(k,k'):=- Y e t (fc)e_ t (£/), 

t6{ —,+} 

i2(fc) := / R(k,k')dk', 


(5.34) 


with t±{k) given by (I5.26p . The remainder terms satisfy 


sup sup V'||r^ ) (A,p)|| I , 2 ( X ) <+cxi (5.35) 

e£(o,i] \ei,\ P \<M j=1 

for any M > 0, compact interval / C (0, +oo). 

From the first equation of the system (I5.32p we get 

D^w e = eW® + | 7 Y *M~ L) +qe, (5.36) 

4G{-, + } 


where 


D ^ \= eX + 2 'yR + ie5 e u}, 


( 5 . 37 ) 
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and q e : = Eti with 


g £ (1) :=e 2 f« 


qf ] ■■= -7 Cu e!+ , qf ] := -ie^R'pu e -. (5.38) 


Computing w e from (15,36j) and then multiplying scalarly both sides of 
the resulting equation by ye M 1 G { —, +} we get the following system 


iw 


W 


37e_e H 
2 D& 


dk — ^—w[ ^ f —dk 
2 e L DU 


= ye 


e,W f 


(o) 


DU 


-dk + 7 


Dh) 


dk , 6 G {—, +}. 


Adding sideways the above equations corresponding to both values of 
l and then dividing both sides of the resulting equation by e we obtain 


a^wi +) 


■a^f 1 (w^U 


where 


w 


<->) = ¥ 


4y f RW f 


(o) 


dk + ^ [ dk , (5.39) 


47 


<K^,p) ■= 

a+(X,p ) := 7 


Dh) ' 3e _ 

' x+i U Rdk 


< T DU 

A + iS e uj 

DU 


(5.40) 


e + dk. 


Equality (15.31 jl is then a consequence of the following. 

Proposition 5.5. For any J G S such that J(y, k ) = J(y) and A > Ao 
we have 


lim 

£^ 0 + 


(f W^Jdpdk - f TJdpdk) = 0, (5.41) 

\J RXT ^ ' JRXT / 

In addition, for any M > 0 and a compact interval I C (0, + 00 ) 

a+ (A, p) (w { e +) (A ,p)~ vjU ] (A, p)) = 0 , (5.42) 


lim sup 

e_>,0 + \ei,\p\<M 


and 


lim sup 

e_>0+ \ei,\p\<M 


a ( w(*,P) 


2A 

T 


0 , 


lim sup 

e_>0+ \ei,\p\<M 


47 

[ R f.dk 

3e 

/ T DU 


0 . 


(5.43) 

(5.44) 


Proof. Equality (15.41 ft follows from (12.27ft and (15.23ft . Formulas (15.42ft 
and (I5.43P can be easily substantiated using the Lebesgue dominated 
convergence theorem. To argue (j5.44j) it suffices only to prove that 


lim sup 

e_>0+ XeI,\p\<M 



Rq, 


(?) 


Dk 


-dk 


0, j — 2 , 3. 


(5.45) 
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Note that 


47 /' Rq, 


(3) 


3e L 


■dk 


< 


47 2 |p|Halloo f R\u e -\ 


| 


dk<C \u e -\dk 


and (15.45ft for j — 3 follows from (I5.28jl . 

In the case j = 2 note that J T dk = 0, therefore 

j£L f ^ [ (A + i^)j 2) ,, 

3e ,/ T DW 3 J T DU 

Since k ha 5 e uj{k,p) is odd for each p we get that the latter integral 
equals 


4 : = 


27 f A(eA + 27 R) + e(8 e ojy 


\DU\ 


Cu ei+ dk. 


(5.46) 


u 


,-i y 

2 ^ 

tG{ —,+} 

47 f Ru e ,+ 


[A(eA + 27 R e ) + e(8 e u) 2 ] dk 


\DU\ 


e ’+./ T \DU\ 2 

[A(eA + 27 R e ) + e(5 e <^) 2 ] dk, 


with 


U eJ(^iP) ( M e,t(A), i,i' e { — ,+}• 


We conclude therefore that there exists a constant C > 0 such that 


|/,l <C 7 e Y, 
(■G{ —, + } 


u 


e ,+1 


MtW) 

\DU\- 


-dk + 76 


R\u e , + \(8 e ujy 

\DU\ 2 


dk. (5.47) 


Denote the terms appearing on the right hand side by li l> and /i 2 \ 
respectively. Since |D^^| 2 > 2e\^R we conclude that there exists C > 0 
such that 


4 (1) < C y |< + 1 and y <C\\u £ , + \\ lHt) . 

tG{—,+} 

Therefore (I5.45jl for j = 2 follows directly from (15.2811 . □ 


5.6. The limit of energy functional at the superdiffusive time 
scale. We have shown in Theorem 15.11 that the evolution of the tem¬ 
perature profile takes place on a scale longer than the hyperbolic one. 
In fact the right time-space scaling is given by (te~ 3 ^ 2 , xe” 1 ) as it has 
been shown in Theorem 3.1 of [5j. 
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Theorem 5.6. Suppose that the distribution of the initial configu¬ 
ration (t, p) satisfies condition (I2.26p . Then, for any test function 
J G C“([0, Too) x R) we have: 


p+oo f t \ r+°° r 

,'TA £ f J(( - ex)E ‘ c *\jn) it = l f 

Here T(t,y ) satisfies the fractional heat equation: 

d t T(t,y) = -c\A y \ 3/4 T(t,y) 


T(t, y)J{t, y)dtdy. 

(5.48) 

(5.49) 


with the initial condition T(0,y) = T(y), given by (j2.27j) and 

[ a "( 0)] 3 / 4 

C 2 9 / 4 (3 7 ) 1 / 2 ' 


(5.50) 


6. Proof of Theorem 13.11 and Corollary 13.21 


Since the dynamics is linear, the solutions of (\2.7\i are of the form 
*x(t) = t' x (t) + t "(f) and p x (t) = p' x (t) + p"(t), t > 0, x G Z, 


where (r'(t),p'(t)) and (r"(f),p"(t)) are the solutions corresponding to 
the initial configurations (t(,,p(,) and (t",p") under the dynamics (12.71) . 
see Section 13.11 The convergence of the v* and p components of the 
vector tti(f) from Corollary 13.21 is therefore a direct consequence of the 
conclusions of Theorems 15.11 and 14.71 

The statements concerning the asymptotics of the thermal and phononic 
components of the energy functional contained in (13.9jl and (13. lOj) also 
follow from the aforementioned theorems. To prove (13.lip , finishing in 
this way also the proof of Corollary 13.21 note that 


^ ) ^th,x 


e J + C ph,. 


- +2t x 


- ,P - 


~) 


t 

( 6 T) 

where e t h,x(t) and e p h, x (^) are defined in (13.41) . Given two conhgurations 
(r( J ), p( J )), j = 1,2 the ’’mixed” energy functional is dehned as 


p<V 7 p«) := 


( 6 . 2 ) 


Here q xx , are computed from (12.91) for the respective conhgurations 
j = 1,2. Formula (13.lip is a simple consequence of the following. 


Lemma 6.1. Under the assumptions of Theorem \3.1\ we have 



(6.3) 
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for any J G Cq°(r) and t > 0. 


Proof. Let (see (12.20[) ) 

(0 := X X (0 - r (^ ex )i 5 ep" Q := p" Q - p(t, ex). 
Using Cauchy-Schwartz inequality we can estimate 



The first factor on the right hand side stays bounded, due to the con¬ 
servation of energy property of the dynamics, while the second one 
vanishes thanks to (14.381) . 

We can conclude therefore that 



0, (6.4) 


for any t > 0, Equality (16.81) would follow, provided we can show that 

f U “V 6 X! J(ex)E ( z x |V Q ,p' Q ; r e(t),p e (t) S j = 0, (6.5) 

for any J G Co°(r) and t > 0, where 

(r e (t),p e (t)) := (r(t,ex),p(t,ex)) xez . 

The latter however is a direct consequence of (15.15j) . 

□ 


7. Proof of Theorem 13.31 


Using the notation from Section [6] we can write the analogue of (16. ip 
at the time scale t/e 3//2 . Thanks to (14.401) for any J G U^°(r) we have 


lim E j (c£) E ee p h,x 

x 



= 0, 


(7.1) 


The respective time-space weak limit of e t h,x (t/e 3//2 ) can be evaluated 
using Theorem 13.31 Finally, to finish the proof of Theorem 13.31 we need 
the following analogue of Lemma 16.11 
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Lemma 7.1. Under the assumptions of Theorem \3.3\ we have 



(7.2) 


for any J G C^°(r) 

Proof. Using Cauchy-Schwartz inequality we can estimate 



The first factor on the right hand side stays bounded, due to the con¬ 
servation of energy property of the dynamics, while the second one 
vanishes thanks to (17.11) . This ends the proof of the lemma. □ 


]. Proof of Theorem 13.41 


From the definition of the normal modes, see (13. 1 71) . we conclude 
that 


df £ \t,k) = \ ± (l-e- 2 “)VEf (±, M) 


( 8 . 1 ) 


D±(l - e 2 ") 2 f (± >((, k) + 0(e 3 (k)) 1 dt + dM,(k), 


where 

dM t (k ) := 2i y 1 ^ 2 J r(k, k')p(t, k — k')B(dt, dk'), (8.2) 

with B(dt.dk) Gaussian white noise in (t,k), satisfying 
B*(dt, dk) = B(dt, —dk), 


E [B(dt, dk)B*(dt ', dk!)] = S(t - t')5(k - k') 
and D± are given by 

37 


2 


(8.3) 
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Let 

(t, k ) := eE e f (± ^ ek^j , k G e _1 T. 

It satisfies 

Lf<±> ((, k) = ±1(1 - e- 2 ")?« (0, k) + 0(e). (8.4) 

An elementary stability theory for solutions of ordinary differential 
equations guarantees that for any T,M > 0 we have 

?W (t, k) = e ±2Mt (l + 0(e)) fW (t, fc), \k\ < M, \t\ < T. (8.5) 


After a straightforward calculation we obtain that the left hand side of 
(I3.20P equals 

lim [ J(p)f ■ ± ’ > (t, —p) dp — lim / (0, — p) dp, 


e->0+ 


e-)-0+ 


which tends to the expression on the right hand side of (13.20H , as 
e —> 0+ for any J G 0£°(r). 

Let 


t 


ef \t,k ) :=fW -,fc exp - 6 (l-e-“) 




, l e {-, +} 


for k G e X T. For any T, M > 0 and t G {—, +} we have 


do 


(0 


dt 


-(t, k) = -^(1 - e M ) 2 d«(U fc) + O(e), |fc| < M, t G [0, T]. (8.6) 


Using again the elementary stability theory of ordinary differential 
equations we conclude that 


lim sup 

e ^ 0+ t£{0,T],\k\<M 


^ ] {t,k) - fW (0, k) e 


47r 2 t.D z ,|/c|' : 


0 , , +}• 


For any J G Cq°(r) we can write that the expression under the limit 
on the left hand side of (13.21)1 equals 

lim^ j J{p ) exp dl'nipy/rl- j — p^j dp 

= lim J j(p) exp j tv ^T (l - e ~ 2Tcitp - 27 Hep) ~ 2 | 6^ (t, -p) dp 

= lim / J(p) exp {2LTjr['Ji 2 p 2 t + O(e)} (t, — p) dp 

e ^°+ Jr 

and the latest limit equals 

[ J(p)f iL \-p)e-^ 2tD ^dp, 

Jr 
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with p u \p) the Fourier transform of p^{y), given by (I3.19p . which 
ends the proof of (13.211) . □ 


9 . Examples 

In the final section we give examples of the initial data that are 
either of thermal or phononic types introduced in Definitions 12.21 and 
[2731 The examples are formulated in terms of the wave function. 


9.1. Non-random initial data. Suppose that (f>(x) is a function that 
belongs to 0 ^°(r). Let a > 0 and let be 5-type measures on i 2 

concentrated at 

;= e ( a ~ l)/2 (f)(6 a x), 16 Z. 


(9.1) 

We have 




Aik) = e (a_1)/2 f 0(p) 

Jr 

^ exp {2t rix(e a p — k)} 

_x£Z 

dp. 



Using Poisson summation formula (see e.g. p. 566 of [6]) 


^exp{ifrr£} = ft f G R ’ ( 9 - 2 ) 

xez ' ' xez ' ' 

(understood in the distribution sense) that holds for any b ^ 0, and 
the fact that 0(fc) is rapidly decaying we conclude 

Uo « (£) ■ 


9.1.1. Case a = 1 - macroscopic initial data. Note that then 


lim sup 

e—>0+ 



e(\m\ 2 ) 

' / Me 


l dp < Too, 


(9.3) 


where <j)(p) is the Fourier transform of (j>{y). The data is of phononic 
type, in the sense of Definition 12.31 with the macroscopic profile given 
by 4>(y). On the other hand, condition (12.261) fails, as can be seen from 
the following computation: 



e(\m\ 2 ) 

\ / n< 


r „ 1 r 1 / 2 

e 2 JjA(k)\ 4 dk « - j ^ 
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9.1.2. Case a G (0,1) - oscillating (but not too fast) data. In this case 
one can easily verify that condition (12.25j) holds, but again condition 
(I2.26P fails. However, the rate of the blow-up is slower than in the case 
of the macroscopic data. Indeed, 


dk 


e(\^(k)( 


Me 


= e 


\^e{k)\ A dk 


, 1 /( 26 “) 


'- 1 /( 26 “) 


\cf)(k) j [ dk ~ — / 1 4>{p) | A dp. 


On the other hand, it is easy to verify that the macroscopic profile for 
the initial data vanishes but 




° = O^i 2 )^ = /1^) \ 2dk - 


Hence the family (/r e ) e> o is neither of phononic, nor thermal type. 


9.1.3. Case a = 0 - microscopically oscillatory data. We have 
V> e (/c) = e -1 / 2 ^ 4>(x) exp {—2irixk} = e -1//2 0(/c), 

xGZ 

where, the periodized Fourier transform of fi{y) is given by 

4>(k) = ^fi{x + k). 

xez 


In this case condition (I2.26P holds, since 


e(\m( 


J dk 

The data is of thermal type. 

9.2. Random initial data. 


Me 


= / \4> (k) \ 4 dk. 


(9.4) 


9.2.1. Modified stationary field. Assume that {j]x)x ez is a zero mean, 
random stationary field such that E1 770 1 2 < + 00 . We suppose that its 
covariance can be written as 


r x = E ( rjxVo ) = / exp {2nikx} R(k)dk, x G Z. 


(9.5) 


where R G C(t) is non-negative. Given a > 0 and fi(x) G Cq^R) 
dehne the wave function as 


:= e*' 0 1 ^ 2 (j)(e a x)rix, iGZ. (9.6) 


One can easily check that condtion (I2.25P holds. We show that both 
micro- and macroscopically varying initial data satisfy condition (I2.26p . 
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For a G (0,1] (the oscillatory case) we have 




c{a-\-V) / 2 


k-i 


m), 


(9.7) 


where fj(d£ ) is the stochastic spectral measure corresponding to {jj x ) xgz . 
Then, thanks to (I9.7j) we get 


dk (eE\^ e \k)\‘ 


—2a 


dk 


k-i 


R(£)di 


= e 


—2a 


dk 


k-i 


k-i' 


R{£)R{£')d£d£' 


}' 

(9.8) 

V 




Changing variables £' := £' j e a and k := k/e a we get that the last 
expression equals 


r /■ i/(26“) /•! 

I J—l/(2e a ) J -] 


■l/(2e°) 

ij(^)d£ 

'T [i-l/(2e“) ^-l/(2e“) 


k - 


^(jfc-0 R(e a £')dkd£‘ 


4 , 




which, as 6 -> 0+, tends to 


£( 0 ) / i?(£)d£ 


■(*) 


dk 


2 


Thus, condition (I2.26P is clearly satished by this family of fields. In the 
case a = 0 (microscopically varying initial data) the condition is also 
valid, as then 

^ e \k) = e~ 1//2 J 4>(k — £)fj(d£), 

with cj)(k) given by (19.4p . As a result 


dk \eE\^ e) {k)\ 2 Y = dk l fi(k-£) R(£)d£ 


< Too. 


9.2.2. Locally stationary initial data. Assume that (t Xje , p x , € )xez > e £ 
(0,1] is a family locally stationary random fields over a probability 
space (fl, J 7 , P). By the above we mean the fields that satisfy the 
following: 

1 ) they are square integrable for each e and there exist C£° func¬ 
tions r, p : R —» R, called the mean profiles satisfying 

<,6 : = = r i ex )i p",6 : = Epi e) = P( ea 0> 
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2) the covariance matrix of the held 

^x,e — ^x,ei Pz,e m ~ P:c,e — Px,e; 

is given by 

E KX+*', e ] = A E [Pa^ePL+s^e] = ^22(ex, x'), 

EK.eP'x+x'^] = ^ 12 (ex, x'), = C , 2 i(ex,x'), x,x' G z, 

where Cij : R x Z AR, i,j = 1,2 are functions that satisfy 


[ Cufax' 

i= 1 m' 


)dx 


< +oo. 


(9.9) 


The energy spectrum of the held (t^ e, Pa: e)xez, cf (13.311 , equals 

a(k) 


txjg(fc) = ^ C 2 2(ex, fc) + 


4s 2 (k) 


^Cn(ex, k) 


where 


Cij(ex,k) = (ex, x') exp {—27rffcx / } , (x, fc) G Z x T, f,j = 1,2. 

x' 

We can write 

e 2 [{tv' e ) 2 dk < Ce 2 EE/ C'jj(exi, k)C*j(ex 2 , k)dk, e G (0,1], 

J T • i _ _ T 


J=1 rri,iT2 ’ 


for some constant C > 0. By an applications of the Plancherel identity 
we conclude that the right hand side approximates (PC*, as e <C 1. 
Thus condition (13.ip holds in this case. 


9.2.3. Local Gibbs measures. Another important example of random 
initial data is furnished by the local Gibbs measure. Given the prohles 
of temperature /3(ex) _1 , momentum p(ex), and tension r(ex), where 
/3 -1 ,p, r G C“(r) and /3 _1 > 0 we dehne a product measure analogous 
to (11.61) . in which the constant prohles are replaced by slowly varying 
functions. These measures are formally written as 

dv x , e ■= ]^[exp j-/3(ex) - p(ex)p x - t(cx)x^\ - £(A(ex)) j dx x dp x , 

(9.10) 

where A(x) = (/3(x),p(x),r(x)) and Q(-) is an appropriate Gibbs po¬ 
tential that normalizes the respective measure. In order to make the 
above ’’definition” rigorous one would have to consider the Gibbs mea¬ 
sure in question as solutions of the respective DLR equations. We shall 
omit that issue by dealing only with local Gibbs measures, i.e. the case 
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when /?(?/) _1 is compactly supported so the relevant measure is defined 
on a finite dimensional space. On the sites where /3(ex ) _1 = 0, we let 
the corresponding exponential factor in (19.10[) be a delta distribution 
concentrated at the point (0,0). The corresponding profile of volume 
stretch r(ex) and temperature f3(ex)~ l are given by analogues of re¬ 
lations (11.9p with r and u appearing there replaced by the respective 
slowly varying functions. 

The natural decomposition in thermal and mechanical initial condi¬ 
tions is now given by 

= r(ex) + x' x , p* = p(ex) + p^ 
where (t^p^x are distributed by 

dv'x, e := “ 0(0(ex), 0,0 )} dx' x dp' x (9.11) 

X 

and t' x is the energy at site x of the configuration (t(,, p' x ) x& z ■ I 11 this 
case we have 

C 2 :2 (ex,x') = S x ',o((p' x ) 2 )v' Xe = o/3~\ex). (9.12) 


9.2.4. Nearest neighbor interactions. Consider first the nearest neigh¬ 
bor case, i.e. when z' x = V{x' x ). Then, 


C*ii(€x 5 x) )*/ 


J r 2 e ~hh x )V(r)^ r 


"A.e J r e -P( ex )V(r)(l r 

Let us assume that there exists c* > 0, for which c*r 2 < V(r). Then, 


0 < Cn(ex, x') < 



f R V(r)e hhx)v{r)^ r 

f R e -/3(ex)V(r)^ r 



Here it (ex) := u(Q,/3(ex)) and u(r, (3) is the internal energy function 
defined in (11.911 . see e.g. (2.1.9) of (2], Condition (19. 9 p is satisfied, once 
we assume that 


/3 ( y)dy < +00 and / u(y)dy < + 00 . 


9.2.5. Gaussian, locally Gibbs measures. We assume that supp/3 " 1 = 
[— K,K] and that fi~ l (y) > 0 for \y\ < K. In addition, we suppose 
that the sequence (a x ), besides satisfying conditions al)-a3) and (12.ip . 
is compactly supported, i.e. there exists a positive integer £ such that 
a x = 0 for \x\ > £ and a x < 0 for 0 < |x| < l (obviously in light of 
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(12.lj) we have cto > 0). In this case, the formal expression (19.lip is a 
probability measure on G( z ) that equals (we omit writing primes) 

V\,e := Pe,K 5 k c , (9.13) 

where §k c is the Borel probability measure on the space of sequences 
(On >e~ 1 K, concentrated on the sequence x x = 0, \x\ > K/e. Mea¬ 
sure [i ( .k is Gaussian on the Euclidean space corresponding to finite 
sequences ( p x )\x\<e~ 1 K whose characteristic functional equals 

6XP I ~\ a x,x'PxPx' > , 

^ \x\,\x'\<e~ 1 K J 

where £ := \cr x ,y\ is the inverse of the symmetric matrix S := [S Xty ] 
corresponding to the quadratic form 

^ J ^ ^ Pi.^')^x—x'^ XlX , j 

\x\,\x'\<e~ 1 K \x\,\x , \<e~ 1 K 

where q X:X > is dehned in (12.9)1 . We claim that there exists c* > 0 
independent of e and such that 

^ ' Sx,x'^x^-x' ^ C* ^ ' G, V (t x ) |a:|<e _1 Jf • (9.14) 

|x|,|x'|<e _1 A |a:|<e _1 A 

According to our assumptions there exists /3* > 0 such that j3(y) > /3* 
for all y E R. We can write therefore 


inf ^ ^ ' S x x 'X x x x i . ^ 1 

\x\,\x'\<e~ 1 K 


\x\<e~ 1 K 


(9.15) 


> — iu f 

~ 4 


£ 


Qt-rr. — 




q.xv : ^ = 1 

|ir|<e — 


> 


A 


inf i _ 5Z <**-*<&<* : = 1 > 

L x,x' X J 


One can easily verify the following identity (cf (12.3)1 ) 

~l^2»x-x'ql,x' = n I y>{s 2 {k))\x(k)\ 2 dk, 


X,X' 


therefore the utmost left hand side (I9.15P can be estimated from below 
by 
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Hence (19.14j) holds with c* as defined above. 

We can write 

= E a x,x' e x(k)e* x t[k), 

|x|,|x'|<e —1 K 

with e x [k) := exp {— 2-Kixk}. Thanks to (19.141) we obtain 

<|e(*)IV. < 4 E KWI 2 = 2</ 't 1) - vfceT - < 9 - 16 ) 

|:c|<e —1 if 

Thanks to (19.12p we conclude that there exists C > 0 such that 

(|p( fc )l 2 )»*,« = E^ i 'Ne = E /rl ( ex ) - p Vee (°> 1 ]- ( 9 - 17 ) 

X X 

Combining (19.161) with (19.17ft we conclude condition (12.261) . 

Appendix A. Proof of Proposition 12.11 
A simple calculation, using the fact that a x = 0, shows that 

E( a * ~ q y ) 2 (A.l) 

x x,y 

for any (q x ) such that r := V*q is square summable. 

Using (1A.1I) , the Plancherel identity and (j2.3[) we can write 

^e x (r,p) = ^ +d(fc)|q(fc)| 2 ^ dk 

= / + Ti<p 2 (s 2 (A:))|r(fc)| 2 ^ dk. (A.2) 

Here r(k), p(k), q (k) are the Fourier transforms of t x , p x , q x respec¬ 
tively. In light of the assumptions made about </?(•) it is clear that the 
utmost right hand side is equivalent to (t x + p x ), so (12.16ft follows. 
Obviously only the lower bound in (12. 1 7ft requires a proof. Note that 

Ew^Etf+Ei^-^- ( A - 3 ) 

X X x,x' 

Here q XjX ' is given by ( 12.91) . By Cauchy-Schwartz inequality for x > x' 
we have 

Qx,x' <(x- x 1 ) E x l" 

x'<x"<x 
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and an analogous inequality holds also for x < x'. Therefore 

^2 I ®-x—x' | t\x,x' — ^ ^ |*£ I |®X — x' \^x" 

x,x' x>x">x' 

+ ^2 \ x - x'\\a x - x :\x 2 x ". 

x'>x">x 

Substituting z := x — x' in the first summation and z x' — x in the 
second we get that the right hand side equals 

J2 z \ a *\ 5Z +^ z \®z\ r *'■ 

z> 0 x>x">x—z z> 0 x+z>x">x 

Denote the first and the second term by I and II respectively. We have 




z >0 


V .// -.// 


x"+ 2 :>a:>a; /y 


< 




, 2>0 


On the other hand 


n = ^2 z \a g J2 1 ^ Z/ 

2>0 x" x+z>x">x \2>0 ) x" 


and, as a result, we get 


+ r2 )’ 


(A.4) 


where C := max {l/2, ^„ >0 z 2 \a z \}. This combined with the already 
shown estimate (12.161) ends the proof of the lower bound in (j2.17[) . 
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